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Extreme events (sometimes also called optical rogue waves), in the form of pulses of extraordinary 
intensity, are easily observed in its chaotic regime if the Fresnel number of the cavity is high. This 
result suggests that the nonlinear interaction among transverse modes is an essential ingredient in 
the formation of extreme events in this type of lasers, but there is no theoretical description of the 
phenomenon yet. We report here a set of experimental results on the regularities of these extreme 
events, to provide a basis for the development of such a description. Among these results, we point 
out here: i) the decay of the correlation across the transversal section of the laser beam, and ii) the 
appearance of extreme events even if the time elapsed since the previous pulse is relatively short 
(in terms of the average inter-pulse separation), what indicates the existence of some unknown 
mechanism of energy storage. We hypothesize that this mechanism is related with the imperfect 
depletion of the gain by some of the transversal modes. We also present evidence in support of this 
hypothesis. 


I. INTRODUCTION 

In recent years there has been a growing interest in 
extreme events (EEs) in various disciplines [T]. The first 
reliable measurements of freak or rogue waves (i.e. un¬ 
expectedly large oceanic waves, as high as 30 m from 
crest to trough) in the beginning of the 90s, were followed 
by intensive research on what had previously been con¬ 
sidered a near-mythological phenomenon. In 2007 Solli 
et al. introduced the concept of optical rogue waves to 
describe large fluctuations in the edge of the spectrum 
of the light generated by the propagation of seed pulses 
in a microstructured optical fibre, thus posing an anal¬ 
ogy between the optical pulses and their oceanic coun¬ 
terparts [2]. The authors did so based on, firstly, the 
L-shaped statistics of the optical events (long tailed dis¬ 
tributions imply the existence of EEs, that, although 
rare, are observed with non-negligible probability), and, 
secondly, the theoretical approach they employed to de¬ 
scribe the results of the experiment. This description is 
based on the Nonlinear Schrodinger Equation (used in 
this case to describe the propagation of pulses in the op¬ 
tical fiber). Rogue waves appear after a modulational 
instability is reached [3]. Since then, similar phenomena 
have been found and studied in a wide variety of op¬ 
tical systems, under the broad denomination of optical 
rogue waves or EEs, meaning, roughly speaking, extreme 
fluctuations in the value of an optical field |1]. How- 
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ever, we must warn the Reader that this does not imply 
a priori a common generation mechanism of these phe¬ 
nomena, or even a common single definition. In fact, 
even in Oceanography, where the term rogue wave was 
originally coined, a unified definition does not yet exist 
. The task becomes more difficult when considering di¬ 
verse physical systems beyond the oceans and analogies 
between these phenomena must be dealt with carefully 
(a thorough discussion on this topic can be found in HD- 
Regardless of their connection with their oceanic coun¬ 
terparts, optical EEs are intrinsically interesting. They 
have been studied in extended optical systems (includ¬ 
ing a linear and a nonlinear experiment -a laser beam 
focused into a perturbed multimode glass fiber, and an 
optical cavity that uses a liquid crystal light-valve as 
nonlinear medium, respectively-) [5], in optically injected 
semiconductor lasers HHg, and in lasers with saturable 
absorbers, both fast (in a Kerr-lens mode locked laser) 
m and slow (in a passively Q-switched all-solid-state 
Nd:YAG-|-Cr:YAG laser) [TT]. For a comprehensive re¬ 
view, we refer the Reader to [12] . In particular, we re¬ 
ported the existence of EEs in a system very similar to 
the latter: a Nd:YV04-|-Gr:YAG laser [T3], a device of a 
wide practical interest. We showed that EEs are observed 
in chaotic regimes with high dimension of embedding, 
high Fresnel number for the cavity and complex spatial 
transverse patterns of the spot. The standard theoreti¬ 
cal approach based on rate equations for a single mode 
[H], though able to describe many of the dynamical fea¬ 
tures of this system, does not predict the existence of 
EEs. This suggests that interaction of transverse modes 
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is a necessary condition for the formation of EEs in this 
system. One of the aims of this contribution is to provide 
further proof for this claim. Another goal is to perform a 
thorough exploration of the regimes in which EEs appear, 
and of the EEs themselves, to guide the theoretical mod¬ 
elling of the problem, as well as to outline key points that 
should be predicted by such model. With these objectives 
in mind, we analyse some of the spatiotemporal features 
of the dynamics of these events. The understanding of 
the mechanism of formation of EEs in this system may 
conceivably lead to its control and, eventually, to useful 
applications. 


In order to quantitatively define EEs, we adopt the cri¬ 
terion that the pulses with a peak intensity exceeding the 
mean by more than 4 times the standard deviation of the 
peak intensity distribution (dcr threshold) are EEs. This 
criterion is somewhat arbitrary; however, we will show 
later that it is appropriate not only because the events 
that fulfill it are among the highest in a particular regime 
of laser operation, but also because they exhibit peculiar 
dynamical features. Besides, we are specifically inter¬ 
ested in those EEs appearing in regimes with long tailed 
histograms, where EEs appear more frequently than in 
Gaussian distributions. In other words, we are interested 
in distributions with a kurtosis higher than 3. 


This paper is organized as follows: in section |TT] we de¬ 
scribe the laser and the methods to record and analyse 
of the time series. Also, the setup to observe transver¬ 
sal coherence domains (IIBl and the spatial correlation 


among different sections of the spot (IICI. In section III 


we discuss the main experimental results, namely: i) the 
plots of intensities in partial sections of the spot show 
that the EEs are not linked to an unique transverse pat¬ 
tern (III A 11; ii) in the dynamical regimes with EEs, the 
two-point time correlation decays to zero at a distance 
nearly half the size of the spot, while in periodic regimes 
such decay is not observed (IIIA2I; iii) the study of the 
time intervals between successive pulses and the return 
maps of peak pulse intensity suggest that the EEs occur 
in a relatively well defined manifold in the phase space 
and, in consequence, that there is some hope to predict 
them (III B I; this claim is supported by the behavior of 
the intensity return maps (III C); iv) heterodyne interfer- 
ograms indicate the existence of domains of transversal 
coherence, confirming the results of (i) (HID); v) a sum¬ 
marizing diagram shows that the regimes with EEs are 
chaotic or hyperchaotic with embedding dimension larger 
than 6, and that they have relatively complex transverse 
patterns (IIIEl. In what follows, we assume the Reader 
to be familiar with the essentials of self Q-switching the¬ 
ory [T5MT7] . 



PM 


Figure 1. Laser setup. LD: pump laser diode,2 W CW @808 
nm; GL: GRIN lens; ND: Nd:YV04 slab (active medium); 
Ml: folding mirror (R = 250 mm); M2: output mirror (plane); 
SA: Gr:YAG crystal, transmission (unbleached) 90%; L: fo¬ 
cusing lens; FP: fast photodiode; GGD: camera for recording 
spot images; PM: power meter, a = 20°; LI = 130 mm; 
L2=75 mm. The position X of the SA is variable, to obtain 
different dynamical regimes. 


II. EXPERIMENTAL SETUP 

A. The laser 

The setup is shown in figure The output of a 2W 
(@808 nm ) CW laser diode is collimated by a GRIN 
lens and focused down to a spot 0.12 mm diameter into 
a Nd:YV04 crystal, 1% doped and appropriately coated, 
mounted on a water- cooled copper heatsink. The V- 
shaped laser cavity has a folding high reflectivity concave 
mirror (R=l m) and a plane output coupler (reflectivity 
= 98%). The mode size varies between the mirrors, with 
the waist near the output coupler. The average output 
power is measured with the power meter PM placed af¬ 
ter the output coupler. A solid-state saturable absorber 
(Cr:YAG crystal, 90% transmission) is placed between 
the folding mirror and the output coupler at a variable 
distance X. By varying its position, the mode size at the 
saturable absorber changes and hence the condition of 
saturation. As the absorber is displaced along the arm 
of the cavity, different dynamical regimes, such as pe¬ 
riodic behavior (periods 2, 4, 6), a period-three stable 
window m and chaotic regimes with and without EEs 
m are observed. 

One of the output beams at the folding mirror is 
focused into a pin fast photodiode (100 ps risetime), 
connected to a PG oscilloscope (PicoScope®6403B: 500 
MHz bandwidth, 5 GS/s, memory of 1 GS). Time series 
of the self-Q-switching pulse intensities with several thou¬ 
sand pulses are recorded. These series are later analysed 
with the TISEAN software package [IHl HO] ? in order to 
calculate the dimension of embedding and the Lyapunov 
exponents. A VGA GGD camera with 60 fps time res¬ 
olution connected to a PG allows to measure the total 
intensity distribution. The whole setup is mounted on 
an optical table with interferometric stability. 
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the spot, while the other (photodiode B) is mounted on 
a micro-metric translation stage, so that it can measure 
the intensity in different regions of the spot. The inte¬ 
grated intensity of the whole spot is at the same time 
registered with the photodiode FP. This scheme allows 
us to measure, simultaneously, the total intensity of the 
spot and the intensity in two specific regions. 


III. RESULTS 


Figure 2. modified Mach-Zehnder interferometer used to 
study transversal coherence of the spot. BSl: beam splitter, 
transmission 96%; MPMl and MPM2: metallic plane mir¬ 
rors; BEl and BE2: 5x beam expanders; BS2: 50/50 beam 
splitter; D: diafragm; S: screen; CCD: CCD camera; B: block. 

B. Interferograms 

In order to study the transversal coherence of the spot 
in different dynamical regimes we use a modified Mach- 
Zehnder interferometer. The laser beam is collimated 
with a convergent lens (F=100 mm) placed at 90 mm 
from the output coupler, and then directed into the in¬ 
terferometer shown in fig. The configuration of the 
Mach-Zehnder is such that it allows us to obtain, by 
means of a 5x beam expander, a superposition of the 
spot with a magnified partial section of itself. This out¬ 
put is in turn magnified (so as to make easier the task of 
discerning interference lines), projected on a screen, and 
registered with a CCD camera. 

C. Measurement of intensities in partial sections of 
the spot 

The comparison of intensities in different partial sec¬ 
tions of the spot (by means of plots of the intensity in 
one section vs. the intensity in the other) is an interest¬ 
ing way to study its transverse dynamics in a pulse to 
pulse evolution, allowing us to determine, for example, 
whether EEs are linked to specific transverse patterns. A 
complementary approach to characterize the transverse 
dynamics is the two point spatial correlation. Unlike the 
plots of intensities in different sections, it provides a mag¬ 
nitude averaged over a whole time series, and it is useful 
to identify the presence of spatial domains and eventually 
detecting spatio- temporal chaos. 

In order to study these features, the loss beam is mag¬ 
nified with a 5x beam expander and then directed to 
a 50/50 beam-splitter. Both reflected and transmitted 
beams are then measured with two photodiodes, which 
are covered by masks with 1mm diameter pinholes to 
limit the measured section, and connected to the PC os¬ 
cilloscope. One of these photodiodes (which we call pho¬ 
todiode A) is fixed, and is kept measuring the center of 


A. Correlation among different regions of the spot 

1. Intensity plots in different sections of the spot 

We study the evolution of the transverse pattern from 
pulse to pulse, by analysing pairs of time series of peak 
intensities in different sections of the spot. We call lAi 
{IBi) the peak intensity measured by photodiode A (B) 
for the ith pulse of the series. Fig. [^(a-d) shows IB^ vs. 
lAi for different distances d for a regime with EEs, and 
fig-i (e) is a histogram of the total peak intensity. As it 
is expected, for d = 0 (i.e., when both photodiodes reg¬ 
ister the same region of the spot), the relation between 
intensities is linear (see fig. (a)); for d = 1 mm (fig. 
1^ (b)), the plot spreads, showing an irregular distribu¬ 
tion; for d = 3 mm (fig. (c)) it becomes practically 
linear again; and finally, for d = 8 mm (fig. i (d)) it 
spreads even more than in the previous case. Four main 
features arise from these plots: i) the transverse pattern 
changes from pulse to pulse in an irregular fashion (a be¬ 
havior like this has been reported for periodic regimes 
of similar systems in |2I] and |22]; ii) some regions of 
the spot show correlation (e.g. fig. [^(c)), while others 
are almost uncorrelated (e.g. figures^ (b) and (d)); iii) 
EEs appear in different zones of the plots: this implies 
that the spatial configuration of the spot changes from 
one EE to another, i.e. the EEs are not associated to 
a single transverse pattern. Moreover, they occur in a 
wide range of values both for lA and IB, i.e.: some EEs 
correspond to a high intensity value in section A and low 
in section B, others to viceversa, and others with a low 
or high intensity intensity in both sections. This means 
that, in the case that EEs were related to specially bril¬ 
liant regions, these would change position from one EE to 
the other, and would not be associated to a specific loca¬ 
tion in the transverse pattern; iv) despite the seemingly 
random distribution in figures (b) and (d), some regu¬ 
larities arise from the plot. E.g., most of the events are 
contained within two approximately straight lines: this 
might be an indication of the predominance of two differ¬ 
ent transverse configurations that appear with a higher 
frequency than others, but with different total intensities 
each time (i.e. every point in one of the lines would be 
associated to the same spatial configuration as the others 
in the same line, but with a different total intensity). 
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Figure 3. Peak pulse intensity measured in the photodiode 
B (IBi) vs peak pulse intensity in photodiode A (lAi) for 
a chaotic regime with EEs, for different values of horizontal 
separation d: (a) d = 0 mm, (b) d = 1 mm, (c) d = 3 mm, 
(d) d = 8 mm. Thick (red online) dots indicate EEs; (e) 
histogram of total pulse intensity for the series, the kurtosis 
of the distribution is 6.02. The data belong to the same time 
series of the regime depicted with blue circles in fig. Thus, 
the corresponding laser spot is the one shown in the left of 
the lower row of that figure. 


2. Spatial correlation 

We calculate the two-point spatial correlation for pairs 
of peak pulse intensity time series (corresponding to dif¬ 
ferent pairs of sections of the spot). It is defined as: 


C{IA, IB) 


N 


_ 

N N 






( 1 ) 


where IA {IB) is the average intensity over the whole 
time series lA (IB). 

Fig. shows the values of C{IA,IB) for two different 
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Figure 4. Up: spatial correlation among different sections 
of the spot as a function of their transversal separation d for 
two different chaotic regimes with EEs (regime 1, depicted by 
circles) and (regime 2, depicted by squares) and a periodic 
regime, with period 6 (diamonds). Down: laser spots for the 
regimes in the upper plot: from left to right, chaotic regime 1, 
chaotic regime 2, periodic; the arrow indicates the direction 
in which the different regions of the spot were sampled, its 
origin shows the position of the fixed photodiode. 


chaotic regimes with EEs which we call regime 1 (cir¬ 
cles) and 2 (squares) and a periodic regime (diamonds). 
The different regimes are obtained by adjusting the po¬ 
sition of the saturable absorber inside the laser cavity, as 
explained in section |II A| The x-axis corresponds to the 
transversal distance d between the sections observed by 
both photodiodes. 

The line depicted by circles is typical of the regimes 
with EEs. The coefficient C{IA,IB) remains close to 1 
for nearly half of the spot size, and then drops to almost 
0 near the edge of the spot. This behavior differs from 
that observed in spatiotemporal chaotic systems with a 
large number (of the order of tens) of modes (which can 
be roughly estimated in optical systems by the Fresnel 
number), characterized by a rapid exponential decay of 
the correlation [23] [23], and/or a narrow peak with a 
short correlation length The behavior observed in 

our laser is consistent with the existence of a relatively 
few transverse modes defining domains of correlation that 
spread through large areas of the spot from pulse to pulse. 

Some of the regimes with EEs show a non monotonous 
behavior, as the one represented by squares in fig. a lo¬ 
cal peak away from the center, before decaying for larger 
d. This does not contradict our previous statement, 
but implies these regimes are dominated by transverse 
patterns different from those involved in the aforemen¬ 
tioned ones, arguably occupying areas with more com¬ 
plex shapes. This is confirmed by the shape of the whole 
spot registered by the CCD camera. 
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For comparison, we also show in fig. [^the spatial cor¬ 
relation for the case of a periodic regime with period 6 
(diamonds). The correlation decays at first, but near the 
edge of the spot it rises again to a value close to 1. In 
other periodic regimes, the spatial correlation remains 
near 1, except in a few narrow valleys. 


B. Temporal intervals between succesive pulses 

Figure]^ shows, for a typical chaotic regime with EEs, 
plots of the intensity of each pulse against the temporal 
interval between that pulse and the previous one (At—, 
upper graphic), and between that pulse and the next one 
(At-|-, lower graphic). The intensity is scaled so that its 
average value is 100 in arbitrary units. Note that both 
plots exhibit a remarkable regularity: the width of the 
range of values of At— and At-I- associated to extreme 
events is far narrower (roughly 2 and 5 /is respectively) 
than that related to, e.g., average pulses (28 and 25 /rs). 
This is evident for the case of At—. This means that: i) if 
one knows that the next pulse is going to be an extreme 
one, then one is able to predict when it is going to occur; 
ii) once an EE has happened, one can predict the time 
that it takes for the build-up of the next pulse (which is 
most probably not an EE) to appear. If an average pulse 
is considered instead, none of these predictions can be 
made. 

Moreover, the plots provide some hints on the mech¬ 
anism of generation of EEs. The first one, given by the 
upper plot, is somewhat counterintuitive: in a simplistic 
approach, one would expect that, the higher the peak 
intensity, the higher the At—, to allow a longer time to 
accumulate energy from the CW. However, the AI—ee 
(i.e. At— intervals related to EEs) are not particularly 
long; in fact, they are closer to the average interval than 
to the maximum At—, since AI—ee — 36 /is, while the 
longest intervals have a value of almost 48 /rs. This sug¬ 
gests that an average pulse does not totally deplete the 
energy available in the gain medium, but instead, that 
it leaves some energy stored in certain regions (possibly 
due to spatial “hole burning”). This stored energy allows 
a pulse to be an EE in spite of the relatively short time 
elapsed since the last one. Therefore, EEs can be thought 
as pulses that are more efficient in extracting the energy 
stored in the gain medium. The lower plot shows that 
the range of values of At+ associated to EEs is among 
the highest of the set, i.e. the time until the next pulse 
after an EE is always among the longest times that can 
be expected for a particular time series (typically At-|- ^ 
40 /is, while the average At+ is 33 /is). This means that 
EEs efficiently deplete the gain, so that the time until 
the next pulse is necessarily long, to allow a new accu¬ 
mulation of gain. 

The plots in fig. [^show what happens in the immedi¬ 
ate vicinity of a given pulse (i.e. one pulse before or one 
after). It is also convenient to study a wider time scale. 
Fig.i shows, for a time series with 112 EEs, a superpo- 
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Figure 5. Plots for peak pulse intensity as a function of 
interpulse time interval for a chaotic time series with EEs. 
Upper (lower) figure shows the peak intensity of each pulse 
as a function of At— (At-|-), i.e. the temporal interval be¬ 
tween said pulse and the previous (next) one. The horizontal 
line indicates the rogue wave limit. Peak intensities are ex¬ 
pressed in arbitrary units scaled so that an average event has 
an intensity of 100. 


sition of time traces centered at each of the EEs (upper 
plot) and, for comparison, a superposition of the same 
amount of average pulses (lower plot). Note that all the 
pulses immediately preceding or following an EE (i.e., if 
the EE is indexed as the Nth event, these pulses are the 
(N-l-l)th and (N-l)th events) occur within a fairly nar¬ 
row temporal window, that slowly blurs for farther pulses 
(the temporal window in which the (N-|-2)th and (N-2)th 
pulses occur is wider than that of the (N-l-l)th and the 
(N-l)th events, and so on). On the contrary, for average 
pulses (lower plot) there is not regular behavior even in 
the immediate vicinity of the pulse. This suggests that 
the trajectory in phase space corresponding to an EE is 
confined to a relatively well defined manifold. This re¬ 
sult gives some hope that EEs can be predict with time 
enough to control them. 
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Time (|as) 


Figure 6. Upper plot: superposition of 112 time traces 
centered at each of the extreme events of a chaotic regime 
of operation (the horizontal line shows the rogue wave limit); 
lower plot: superposition of time traces centered at 112 near- 
to-average events for the same regime of operation as in the 
upper plot. 


C. Return intensity maps 


Return maps, or plots of the intensity of the peak in¬ 
tensity of the (n-l-l)th pulse vs that of the nth pulse, pro¬ 
vide useful insight. This is shown in fig. [^for the regime 
1 of fig. in Once again, the EEs show a specific and reg¬ 
ular behavior: they appear related to two relatively well 
localized “tongues” emerging from a central bunch at a 
value of ~ 80. They are similar to those corresponding to 
the plot of intensity vs interpulse time interval in fig. 
Therefore, EEs are preceded and followed by relatively 
low intensity pulses (slightly below the average), and the 
intensity of those pulses is restricted to a narrow range of 
values. This means, once again, that the evolution of the 
peak intensity of the pulses surrounding an EE is quite 
repetitive. The fact that the pulse preceding the extreme 
event is close to average supports the picture described 
in section UlIBI 

In summary, according to the results presented in sec¬ 
tions |llT^ and |III^ the EEs are not merely high inten¬ 
sity pulses. They seem to follow a dynamics that is differ¬ 
ent from that of the average Q- switch pulse in the tem- 
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Figure 7. Return intensity map for a chaotic regime of 
operation with extreme events. The vertical and horizontal 
line indicate the rogue wave limit. 


poral series. Furthermore, the dynamics of EEs seems to 
be more regular than that of the average pulses. This 
enforces the idea that there is a deterministic mechanism 
behind the formation of EEs. 


D. Interferograms 


Fig. HI shows the laser spot and heterodyne interfer¬ 
ograms, obtained as explained in section IIB for two 
different dynamical regimes: a chaotic regime without 
rogue waves (upper row), and a chaotic regime with EEs 
(lower row). In the former regime, fringe patterns are 
observed all across the spot. In the latter regime, in¬ 
stead, fringes are observed in the region that is expanded 
on the reference beam and its immediate neighborhood. 
Outside this region, fringes blur and disappear, indicat¬ 
ing that the coherence is lost. This result suggests the 
existence of coherence domains and provides further sup¬ 
port to the transverse mode interaction hypothesis. It 
is important to remind that, because the transverse pat¬ 
tern varies from pulse to pulse (as it was shown in section 


III A 11, and due to the long exposure time of the CCD, 


neither the spot nor the interferogram correspond to a 
single pulse transverse pattern, but rather to a superpo¬ 
sition of many different ones. We foresee repeating this 
experience using an ultrafast camera in order to distin¬ 
guish pulse to pulse spots and interferograms. 
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Figure 8. Each row corresponds to a different regime; up¬ 
per row: chaotic regime without extreme events; lower row: 
chaotic regime with extreme events; left column: laser spot; 
right column: interferogram. The region marked with a black 
line in the lower right figure indicates the region where the 
interference fringes blur 
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Figure 9. Representation of all the recorded dynamical 
regimes in this study having a measurable embedding di¬ 
mension in terms of the number of lobes of their spot (hor¬ 
izontal axis) and their embedding dimension (vertical axis). 
The color and shape code indicates whether the series are 
not chaotic, i.e. periodic (full black squares), chaotic without 
extreme events (hollow blue circles), hyperchaotic without ex¬ 
treme events (full blue circles), chaotic with extreme events 
(hollow red stars) or hyperchaotic with extreme events. 


IV. SUMMARY 


E. Dynamical and spatial complexity 

Fig. ^ summarizes the features of all the time series 
recorded in this study having a measurable embedding 
dimension. They are placed according to their spatial 
complexity, which is quantified by the number of “lobes” 
in the spot and their value of embedding dimension. A 
color and shape code indicates whether a particular se¬ 
ries is non-chaotic, chaotic (one positive Lyapunov expo¬ 
nent), hyperchaotic (more than one positive Lyapunov 
exponent), and whether it displays EEs or not. It can 
be seen, as a general behavior, that chaotic and hyper¬ 
chaotic regimes have a “large” dimension of embedding. 
Here, “large” means higher than the value of 4 predicted 
by the standard theoretical model based on rate equa¬ 
tions for a single mode na. The time series with EEs 
tend to display a large number of lobes, too. This im¬ 
plies that regimes with EEs are associated with a dynam¬ 
ical behavior which is complex both in space and time. 
On the other hand, chaotic with no EEs and periodic 
regimes tend to concentrate in the lower left quadrant 
of the graph, which implies they have a lower degree of 
complexity (lower dimension of embedding and spatially 
simpler spots). This figure enlarges the results presented 
in [IH]. 


In this paper, we report a series of observations aimed 
to guide the way to a (still missing) theoretical expla¬ 
nation of the formation of EEs in all-solid-state, self-Q- 
switched lasers. 

In the regime with EEs, the spatial correlation de¬ 
cays to zero in a distance comparable with the spot size. 
This is not observed in periodic (therefore without EEs) 
regimes. No decay associated with a narrow peak is 
observed as it was, instead, in other extended spatio- 
temporal systems. This is consistent with the idea that 
relatively few transverse modes are involved in the forma¬ 
tion of the EEs. The plots of the intensities recorded in 
the two detectors, or IB vs lA, show that the transverse 
patterns change in an irregular way from pulse to pulse, 
and that the EEs are not linked to a specific transverse 
pattern. 

From the analysis of the interpulse time intervals we 
see that: i) the values of At—EE {At+Es) are typically 
contained in a narrow temporal span; ii) AI—ee are next 
to average, contrarily to what might be expected; iii) 
At+EE, on the other hand, are among the longest in¬ 
terpulse intervals. The result (ii) is consistent with a 
scenario where typical pulses do not totally deplete the 
accumulated gain in the active medium, but leave some 
energy stored, so that the following pulse does not need 
a build-up time proportional to its intensity. This energy 
storage is presumably related with spatial hole burning 
in the active medium. The result (iii) suggests that the 
deep cause of the EEs is simply that they are particularly 












efficient in depleting the gain, through a mechanism that 
is to be elucidated. All these results, together with the 
regular behavior exhibited in fig. suggest that there 
exists a deterministic mechanism in the formation of the 
EEs, and therefore that there is some possibility of pre¬ 
dicting and controlling them 

The interferograms have shown the existence of do¬ 
mains of coherence compatible with a dynamics ruled by 
modes interaction. 

Finally, the summarizing diagram Fig. shows that 
the EEs are prone to arise in chaotic dynamic regimes 
with large dimension of embedding and with a spot with 
a complex transverse structure. This result also supports 
the few- mode interaction hypothesis as the basis of the 
mechanism of formation of EEs in this type of lasers. 

We foresee using an ultrafast camera in order to record 
series of individual pulse spots (as well as single spot 


interferograms) and identify the ones corresponding to 
EEs. Another planned course of action is to replace the 
pump laser diode with a VCSEL, which provides a much 
spatially uniform pump mode, and study its effect on the 
formation of EEs. Our ultimate goal is the construction 
of a theoretical model able to predict the dynamics of 
EEs in this system. 

ACKNOWLEDGEMENTS 

This work was supported by the grant FA9550-13-1- 
0120, “Nonlinear dynamics of self-pulsing all-solid-state 
laser” of the AFOSR (USA), the contract PIP2011- 
077 “Desarrollo de laseres solidos bombeados por dio- 
dos y de algunas de sus aplicaciones” of the CONICET 
(Argentina), and the project OPTIROC of the ANR 
(France). 


[1] S Albeverio, V Jentsch, and H Kantz, editors. Extreme 
Events in Nature and Society. The Frontiers Collection. 
Springer Berlin Heidelberg, Berlin, Heidelberg, 2006. 

[2] D R Solli, C Ropers, P Koonath, and B Jalali. Opti¬ 
cal rogue waves. Nature, 450(7172):1054-1057, December 
2007. 

[3] E N Pelinovski and C Kharif. Extreme Ocean Waves. 
Springer Science & Business Media, June 2008. 

[4] J M Dudley, F Dias, M Erkintalo, and G Genty. Instabili¬ 
ties, breathers and rogue waves in optics. Nat. Photonics, 
8(10):755-764, October 2014. 

[5] V Ruban, Y Kodama, M Ruderman, J Dudley, 
R Grimshaw, P V E McClintock, M Onorato, C Kharif, 
E Pelinovsky, T Soomere, G Lindgren, N Akhmediev, 
A Slunyaev, D Solli, C Ropers, B Jalali, F Dias, and 
A Osborne. Rogue waves - towards a unifying con¬ 
cept?: Discussions and debates. Eur. Phys. J. Spec. Top., 
185(1):5-15, July 2010. 

[6] S Residori, U Bortolozzo, A Montina, F Lenzini, and 
F T Arecchi. Rogue waves in spatially extended optical 
systems. Pluct. Noise Lett., 11(01):1240014, March 2012. 

[7] C Bonatto, M Feyereisen, S Barland, M Giudici, C Ma- 
soller, J Rios Leite, and J Tredicce. Deterministic Optical 
Rogue Waves. Phys. Rev. Lett., 107(5):053901, July 2011. 

[8] J Zamora-Munt, B Garbin, S Barland, M Giudici, J R 
Rios Leite, C Masoller, and J Tredicce. Rogue waves 
in optically injected lasers: Origin, predictability, and 
suppression. Phys. Rev. A, 87(3):035802, March 2013. 

[9] J Reinoso, J Zamora-Munt, and C Masoller. Extreme 
intensity pulses in a semiconductor laser with a short 
external cavity. Phys. Rev. E, 87(6):062913, June 2013. 

[10] M G Kovalsky, A A Hnilo, and J R Tredicce. Extreme 
events in the Ti:sapphire laser. Opt. Lett., 36(22):4449- 
4451, November 2011. 

[11] A A Hnilo, M G Kovalsky, and J R Tredicce. Extreme 
value events in self pulsing lasers. In 2011 Eifth Rio De 
La Plata Workshop on Laser Dynamics and Nonlinear 
Photonics, pages 1-3, 2011. 

[12] N Akhmediev, J M Dudley, D R Solli, and S K Turitsyn. 
Recent progress in investigating optical rogue waves. J. 


Opt., 15(6):060201, June 2013. The full text of this article 
is available in the PDF provided. 

[13] G Bonazzola, A Hnilo, M Kovalsky, and J R Tredicce. 
Optical rogue waves in an all-solid-state laser with a sat¬ 
urable absorber: importance of the spatial effects. J. 
Opt., 15(6):064004, June 2013. 

[14] D Y Tang, S P Ng, L J Qin, and X L Meng. Determin¬ 
istic chaos in a diode-pumped NdYAG laser passively Q 
switched by a Gr 4-1-: YAG crystal. Opt. Lett., 28(5):325- 
327, March 2003. 

[15] A E Siegman. Lasers. University Science Books, January 
1986. 

[16] W Koechner. Solid-State Laser Engineering. Springer, 
New York, NY, 6th, rev. and updated ed. 2006 edition 
edition, April 2006. 

[17] A Yariv. Quantum Electronics. John Wiley & Sons, 
Incorporated, January 1975. 

[18] G Bonazzola, A A Hnilo, M G Kovalsky, and J Tredicce. 
Optical rogue waves in the solid-state laser with a sat¬ 
urable absorber. In Proceedings of the IEEE, pages 1-3, 
Montevideo, Uruguay, December 2013. IEEE. 

[19] T Hegger, H Kantz, and T Schreiber. Practical im¬ 
plementation of nonlinear time series methods: The 
TISEAN package. CHAOS, 9:413, 1999. 

[20] T Schreiber and A Schmitz. Surrogate time series. Phys- 
ica D, 142:346, 2000. 

[21] M-D Wei. Modeling a Passively Q-Switched Nd:YV04 
Laser with an Application to Spatio-Temporal Instabil¬ 
ity. Jpn. J. Appl. Phys., 49(7R):072701, July 2010. 

[22] M-D Wei, G-H Ghen, and K-G Tu. Spatial and tempo¬ 
ral instabilities in a passively Q-switched Nd:YAG laser 
with a Gr4-|-:YAG saturable absorber. Opt. Express, 
12(17):3972-3980, 2004. 

[23] F Arecchi, G Giacomelli, P Ramazza, and S Residori. Ex¬ 
perimental evidence of chaotic itinerancy and spatiotem- 
poral chaos in optics. Phys. Rev. Lett, 65(20):2531-2534, 
November 1990. 

[24] P G Hohenberg and B I Shraiman. Ghaotic behavior of 
an extended system. Physica D, 37(1-3): 109-115, July 
1989. 



9 


[25] G Huyet, M C Martinoni, J R Tredicce, and S Rica. Number. Phys. Rev. Lett., 75(22):4027-4030, November 

Spatiotemporal Dynamics of Lasers with a Large Fresnel 1995. 



